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Abstract 



The cross section for dijet photoproduction at high transverse energies is presented 
as a function of the transverse energies and the pseudorapidities of the jets. The 
measurement is performed using a sample of ep-interactions corresponding to an 
integrated luminosity of 6.3 pb _1 , recorded by the ZEUS detector. Jets are defined 
by applying a fey-clustering algorithm to the hadrons observed in the final state. The 
measured cross sections are compared to next-to-leading order QCD calculations. 
In a kinematic regime where theoretical uncertainties are expected to be small, the 
measured cross sections are higher than these calculations. 
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1 Introduction 



In photoproduction at HERA a quasi real photon, emitted from the incoming positron, 
collides with the incoming proton. In leading order quantum chromo dynamics (QCD), 
two processes contribute to the photoproduction of jets: the direct process, in which 
the photon couples directly to a parton in the proton, and the resolved process, in 
which the photon acts as a source of partons, one of which scatters from a parton in 
the proton. Beyond the leading order in QCD, direct and resolved processes are not 
distinctly separable. 

The cross section for jet photoproduction is sensitive to the partonic structures of both 
the proton and the photon. In the kinematic regime of the measurement presented in 
this paper, the fractional momentum x at which partons inside the proton are probed 
lies predominantly in the region between 10 -2 and 1CP 1 . At these x values the parton 
densities in the proton are strongly constrained by measurements of the structure func- 
tion F% in deep inelastic lepton-proton scattering The fractional momentum x 7 at 
which partons in the photon are probed lies between 0.1 and 1. For x 7 values above 
0.5 the quark densities in the photon are not strongly constrained by i 7 ^ data obtained 
from 77* scattering at e + e~ experiments Q. 

The investigation presented in this paper aims to constrain more tightly the parton 
densities in the photon at high x 7 , where the contribution from quarks dominates, by 
exploiting their influence on the dijet photoproduction cross section. For this purpose 
the dijet cross section is measured in a kinematic regime where next-to-leading order 
(NLO) QCD calculations are expected to describe the data. It should be noted here 
that jet measurements at the Tevatron Q, although generally in good agreement with 
NLO-QCD, show discrepancies in the comparison of the 630 GeV and 1800 GeV data 
sets. These may be connected to non-perturbative effects, such as a possible underlying 
event 0|. A number of these effects, which may also be of relevance to the present study, 
have been investigated in this paper. 

This paper builds on the improved understanding of jet photoproduction and of compar- 



isons to NLO-QCD calculations gained in previous analyses || - [10 and on a significant 



theoretical effort in the recent past (jlTJ - [24|. Events with two or more high-transverse- 
energy jets are used, one of which is required to have transverse energy greater than 
14 GeV and the second one greater than 11 GeV. A previous jet photoproduction anal- 
ysis Q has shown that for jets with transverse energy greater than 11 GeV, the dijet 
cross section agrees with NLO-QCD predictions, within the experimental uncertainties 
of that analysis. 



2 Experimental setup 

The data used in this paper were collected in 1995 with the ZEUS detector at HERA, 
colliding positrons at an energy of E e = 27.5 GeV with protons at an energy of E p = 
820 GeV, yielding a total CM energy of </s = ^J^E e E v « 300 GeV. The data sample 
corresponds an integrated luminosity of 6.3 pb _1 . 

The ZEUS detector is described in detail elsewhere p5| . The main components used 
in this analysis are the uranium-scintillator calorimeter (CAL) and the central tracking 
detector (CTD). The CAL pRpfl covers 99.9% of the total solid angle and is subdivided 



into forward, barrel and rear parts, covering the pseudorapidity regions 4.3 > r\ > 1.1, 



1 



1.1 > K] > —0.75 and —0.75 > i] > —3.8, respectively^. Test beam measurements 
yield energy resolutions of o(E)/E = 18%/^/E(GeV) for electrons and a(E)/E = 



35%/V#(GeV) for hadrons |7). The CTD |2§ is a cylindrical drift chamber, situated 
in a 1.43 T solenoidal magnetic field, covering the polar angular region 15° < 9 < 
164°. The transverse momentum resolution for full-length tracks can be parametrised 
as ct(pt)/pt = 0.0058^t© 0.0065 00.0014/pt, with p T in GeV. The luminosity collected 
by ZEUS is measured from the rate of the Bremsstrahlung process e + p — > e + pj. A 
three-level trigger system is used to select events online |8,25|. 



3 Definition of the cross section 

The relevant variables for the dijet cross section measurement presented in this paper 
are the following: 

• the transverse energy, E^, the azimuthal angle, (f>> et , and the pseudorapidity, rf et , 
of the jets; 

• the scaled energy transfer from the positron to the proton in the proton's rest 
frame, defined as: 

y = l^-,0<y<l, (1) 
k ■ p 

where q, k and p are the four-momenta of the exchanged photon, the incoming 
positron and the incoming proton, respectively. Neglecting mass terms, y is related 
to the centre-of-mass energy in the photon-proton system, W 7P = y/ys. In the 
photoproduction regime, where the exchanged photon is almost real, y is equivalent 
to the fractional energy of the incoming positron carried by the photon; 

• the fractional longitudinal momentum of the photon participating in the produc- 
tion of the two highest-transverse-energy jets, defined as ||: 



X " ~ 2yE, ' ( ) 



where E^ 2 and rf^ are the transverse energies and the pseudorapidities of the 
two highest-transverse-energy jets; 

• the virtuality of the exchanged photon: 

Q 2 = -i 2 ■ (3) 

The cross section presented in this paper is compared to NLO-QCD predictions. It is 
restricted to a specific set of conditions, to minimise theoretical uncertainties. 

• An asymmetric cut is applied on the transverse energy of the two highest-transverse- 
energy jets. The application of a symmetric cut poses a stability problem for some 



of the available NLO-QCD calculations [12, 13]. 



1 The ZEUS coordinate system is denned as right-handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X axis horizontal, pointing towards the centre of 
HERA. The pseudorapidity is defined as r\ — — ln(tan |), where the polar angle 6 is taken with respect 
to the proton beam direction. 



2 



Symmetrisation of the cross section with respect to the pseudorapidity of the two 
highest-transverse-energy jets has been claimed to remove infrared instabilities 



in the NLO-QCD calculations [14]. This entails analysing each event twice, as 
explained below. 

Jets are defined using the longitudinally invariant ^-clustering algorithm pa] in 



the inclusive mode [16], where the parameter R is chosen equal to 1. This algorithm 
provides a jet reconstruction that is suitable for comparisons between data and 
theory 0. 

The dijet photoproduction cross section presented in this paper refers to events in which 
at least two jets, as defined by the ^-clustering algorithm, are found in the hadronic final 
state. These jets are required to have pseudorapidities between —1 and 2, transverse 
energy of the highest-transverse-energy jet, Ej^j eading , greater than 14 GeV and the 

transverse energy of the second-highest-transverse-energy jet, E^ eamd , greater than 11 
GeV. The cross section is given in the kinematic region defined by: Q 2 < 1 GeV 2 and 
0.20 < y < 0.85. 

This cross section is measured as a function of three variables: E^ leading , rf^ 1 and 
rf^ 1 ■ The cross section is symmetrised with respect to the pseudorapidities of the two 
jets. Every event contributes twice to the cross section, once with rj{ et = rfi^ding ana - 

jet jet i i i • -.1 jet jet j jet jet 

m = ^second and a SeCOnd tlme Wlth m = %econd and V2 = Vleading- 

The cross section is determined for the full range of x° k values and for a direct- 
photoproduction-enriched region with x° bs > 0.75. The cross section as a function of 
the pseudorapidity of the jets is also measured in a narrower band of y values between 
0.50 and 0.85, where the sensitivity to the photon structure is expected to be higher, as 
will be explained in section W. 



4 Comparisons to NLO-QCD 

The measured cross sections are compared to NLO-QCD calculations by four different 
groups: P. Aurenche et al. |l8|], S. Frixione et al. [|^,|l^|, B. Harris et al. p0| a nd M. 



Klasen et al. |21J. These calculations differ in the handling of divergences [22,23]. 
All calculations use the CTEQ4M J2j| parameterisation of the parton densities in the 
proton. The value of Aqcd is chosen to match that of this set of parton distribution 
functions. For the parton densities in the photon three parameterisations are used, 
GRV-HO M,M, GS96-HO pi and AFG-HO §3|. 



In all calculations the renormalisation and factorisation scales are chosen equal to the 
transverse energy of the highest-transverse-energy jet. The variation in the NLO-QCD 
calculations of the presented cross section has been found to be less than 15%, when the 
scales are varied between half and twice this value. 

The NLO-QCD calculations do not include fragmentation. Jets are defined on the ba- 
sis of the outgoing partons. While the momenta of jets at high transverse energies are 
expected to correspond closely to the momenta of the partons produced in the hard sub- 
process, the measured jet cross sections are affected at some level by the fragmentation. 
In a study using the HERWIG 5.9 and the PYTHIA 5.7 Monte Carlo photoproduction 
models, the dijet cross section for jets of hadrons was compared to that for partons pro- 



3 



duced in the two-to-two hard subprocess and in the parton showers which were grouped 
into "parton jets" using the fc^-clustering algorithm. 

In HERWIG the change in the cross section due to the fragmentation was found to be 
less than 10% in most of the present kinematic region. However for events in which one 
jet has rf et < the cross section is reduced by more than 10% due to fragmentation 
and when both jets have rf et < the cross section is reduced by ~ 40%. In PYTHIA 
the reduction of the cross section due to fragmentation is much smaller, but shows the 
same trend. In a related study, presented in reference HERWIG 5.9 was used to 
compare the cross section for jets of hadrons to that for jets of partons, produced in 
the two-to-two hard subprocess. The relative difference between these cross sections 
was found to be less than 20%, except again for events with backward jets where the 
reduction of the cross section due to fragmentation exceeds 20% and is again largest 
(~ 50%) when both jets have r^ 6 * < 0. 

Since the effect of fragmentation on the cross section depends on the Monte Carlo 
model, no attempt was made to correct the data for these effects. Instead, the effect of 
fragmentation is considered as a theoretical uncertainty. 



5 Energy corrections 

Kinematic variables are reconstructed using a combination of track and calorimeter 
information that optimises the resolution of reconstructed kinematic variables |34|] . The 
selected calorimeter clusters and tracks are referred to as Energy Flow Objects (EFOs). 
The use of track information reduces the sensitivity to energy losses in inactive material 
in front of the CAL. However, the energy of particles for which no track information is 
available (e.g. because the energy is deposited by a neutral particle), must be measured 
using CAL information. These energies have to be corrected for the energy losses in the 
inactive material. The conservation of energy and momentum in neutral-current deep 



inelastic scattering events is exploited to determine the required energy corrections [35] 
by balancing the scattered positron with the hadronic final state. This is done for 
data and Monte Carlo event samples independently. The EFOs thus corrected are used 
both for the reconstruction of jets and to determine kinematic variables. Comparisons 
between data and Monte Carlo of kinematic variables, reconstructed using corrected 
EFOs, lead to the assignment of a 3% correlated systematic uncertainty and a 2% 
uncorrelated systematic uncertainty in the transverse jet energies and in the hadronic 



variables |35 |. 



6 Event selection 

After applying the energy corrections described in section ||, dijet events are selected 
from those events triggered by the dijet trigger [|| using the following procedures and 
cuts designed to remove sources of background: 

• The ^-clustering algorithm, in the inclusive mode with R = 1, is applied to the 
corrected EFOs. Events are selected in which at least two jets are found with: 
-1 < j« < 2, E^ leadmg > 14 GeV and E^ second > 11 GeV. 
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To remove background due to proton beam-gas interactions and cosmic showers, 
a cut is made on the longitudinal position of the reconstructed interaction vertex 

- 40 cm < Z vertex < 40 cm . (4) 

To remove background due to charged-current deep inelastic scattering events, a 
cut is made on the relative missing transverse momentum: 

-p= < 1.5\/GeV, (5) 

where Pt and Et are the transverse momentum and the transverse energy of the 
event, calculated on the basis of corrected EFOs. 

The rejection of neutral-current deep inelastic scattering (NC-DIS) events is based 
on the variable y. If a scattered positron candidate with energy greater than 5 GeV 
is found in the calorimeter, y can be calculated from the energy E' e and the polar 
angle 9' e of this positron candidate using the formula: y e i ec = 1 — JIT ~~ cos ^e)- 
These events are rejected when: 

Velec< 0-7. (6) 

The variable y can also be reconstructed from the observed hadronic final state 
using the Jacquet-Blondel formula [36]: yjB = J2{E ~ Pz)/2E e , where the sum 



runs over all corrected EFOs. For all events it is required that 

0.20 < y JB < 0.85 . (7) 

This cut removes unidentified NC-DIS events, for which yjs peaks at 1, and proton 
beam-gas interactions, which mostly have low yjs values. The cuts on y e i ec and 
yjB effectively restrict the range of the virtuality of the exchanged photon to 
Q 2 < 1 GeV 2 , with a median of about 10~ 3 GeV 2 . 

After the application of the described selection criteria, a sample of 8690 events re- 
main. The contamination of this sample due to background processes was found to be 
negligible. 

7 Event characteristics 

Photoproduction events, generated using Monte Carlo programs, are used for the de- 
termination of acceptance and migration corrections and for the study of systematic 
uncertainties. These events are passed through a full simulation of the ZEUS detector 
and undergo the same energy-correction procedure as the data. 

Two leading-order Monte Carlo programs were used to generate dijet photoproduction 



events, HERWIG 5.9 |7|,||] and PYTHIA 5.7 US. Both models use leading-order 
matrix elements, but they differ in the treatment of parton showers, hadronisation and 
the virtuality spectrum of the exchanged photon. No additional process that would 
produce soft or hard underlying events is included in the simulations. 
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Direct and resolved event samples are generated separately. The parton density functions 
used to generate both Monte Carlo samples are CTEQ3-LO [41] for the proton and 
GRV-LO mm for the photon. 

As the Monte Carlo models do not include higher order matrix elements, they are not 
expected to describe the absolute normalisation of the cross section. To obtain the best 
agreement between data and Monte Carlo, the normalisation of the direct and resolved 
contributions is determined from a fit to the measured spectrum. As a result the 
direct contribution of the HERWIG Monte Carlo is scaled by a factor 1.83 and the 
resolved contribution by a factor 1.72. For PYTHIA the direct contribution is scaled by 
1.28 and the resolved contribution by 1.27. When these factors are applied, both Monte 
Carlo models are found to give a reasonable description of various distributions, such as 
the rf et and E 3 ^ spectra. 

The x° bs spectrum for the selected sample of the 8690 dijet events is shown in Fig. |l|, 
where x° bs is determined on the basis of corrected EFOs. The data show a clear peak 
near x° bs ~ 1, attributed, at leading order to a predominance of "direct" events, and 
a tail towards low x° bs values, attributed to "resolved" events. The data are compared 
to the HERWIG and PYTHIA Monte Carlo predictions, including the normalisation 
factors given above. The direct component of the HERWIG Monte Carlo is shown sepa- 
rately. The direct photoproduction events peak at high x° bs values. Therefore, selecting 
events with x° bs > 0.75 yields a sample strongly enriched with direct photoproduction 
events. After application of the normalisation factors described above, the Monte Carlo 
predictions are in good agreement with the data. The shape of the peak at high x^ s is 
best described by the HERWIG Monte Carlo. Given the agreement in this distribution 
and in distributions like the rf et and E 3 ^ 1 spectra (not shown here), the HERWIG Monte 
Carlo sample is used to determine acceptance and migration corrections and to study 
systematic uncertainties. 

In Fig. |2] the transverse energy flow around jets is shown as a function of the distance in 
pseudorapidity At], with respect to the jet axis, integrated over Acfi, between cj>> et — 1 and 
</> je * + l, where cj) is measured in radians. The transverse energy flows are shown in bins of 
E 3 ^ and i*. The jets are strongly collimated, with relatively little transverse energy 
away from the jets. Comparison to the HERWIG predictions shows generally good 
agreement. Only at low x° bs values is the energy flow outside the jet underestimated by 
the Monte Carlo model. Jets in the Monte Carlo are also found to be slightly narrower 
than jets in the data. 

In the kinematic regime of the present analysis, Monte Carlo models that do not include 
a simulation of underlying events are able to describe the x^ s distribution and the 
transverse energy flows, the only exception being the transverse energy flows in the 
lowest x° bs bins, where a small discrepancy is observed. As these distributions are 
considered to be particularly sensitive to underlying events, this indicates that such 
processes play no role in the present kinematic regime. This result is different from 
what was observed in previous photoproduction analyses |^|, in which jets at lower 
transverse energies were studied, where it was shown that the description of the data 
is improved when a model simulating soft or hard underlying events is included in the 
simulations. 
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8 Unfolding and Systematics 



The unfolding of the cross section is done by multiplying the number of events recon- 
structed in each bin by a correction factor determined from the HERWIG Monte Carlo 
sample. This correction factor is defined as the ratio of the number of events generated in 
the bin, Nt rU e, over N rec , the number of events reconstructed in the bin. The systematic 
uncertainty related to the choice of Monte Carlo model for the unfolding is estimated by 
using a different Monte Carlo generator, PYTHIA, to determine the correction factors. 
The HERWIG and PYTHIA Monte Carlo models differ in the treatment of the gener- 
ation of the photon spectrum, the parton showers and the simulation of hadronisation 
effects. Nevertheless, both Monte Carlo models give a reasonable description of the 
data. 

To determine systematic uncertainties in the measured cross sections, several variations 
in the event selection have been studied. The uncertainty in the cross section due to the 
energy-scale uncertainty is estimated by raising and lowering all energies in the Monte 
Carlo simulation by 3% simultaneously. In addition y and the transverse jet energies 
are varied by ±2% separately. 

The systematic uncertainty related to the energy measurement is correlated from bin 
to bin. In the cross section figures presented in the next section, this uncertainty is 
shown separately. All other positive (negative) contributions to the cross section, from 
systematic uncertainties, are added in quadrature to yield the total positive (negative) 
systematic uncertainty. 

9 Results 

The dijet photoproduction cross section is presented as a function of three variables: 
Lading' the transverse energy of the leading jet, and rf^ 1 and t^ 6 *, the pseudorapidities 
of the two jets. Statistical and systematic uncertainties, added in quadrature, are shown 
as thin error bars. Statistical uncertainties alone are shown as thick error bars and the 
uncertainty due to the energy scale is shown as a shaded band. 

9.1 Cross sections for 134 < W 1P < 277 GeV 

The dijet cross section as a function of the transverse energy of the leading jet is pre- 
sented for six different ranges in jet pseudorapidity. These cross sections have been 
determined both for the full x° bs range and for x° bs > 0.75. Numerical values for the 
cross sections and the uncertainties are given in Tab. [I] and |2| The results are plotted 
in Fig. |3] and |j. The dijet cross section falls rapidly with increasing transverse energy of 
the leading jet. The steepest slopes occur when both jets are in the most backward pseu- 
dorapidity bin, — 1 < rfi2 < 0- High x° bs events dominate the cross section at backward 
angles of the jets and at high transverse energies of the jets. This behaviour is expected 
on kinematic grounds, since high x° bs values give access to the highest transverse jet 
energies and to the most backward pseudorapidities. 

The data are compared to NLO-QCD calculations (see section ||). Since the calculations 
from different groups are very similar, as will be shown in Fig. || and |(| only one set of 
calculations is shown here. This set corresponds to the GRV-HO |30| , [3ll| parameterisation 
of the photon structure, which gives the highest cross section. In general, the slopes and 



7 



the absolute cross section are well described by the NLO-QCD calculations. However 
for events with forward jets, 1 < rf-^ < 2, and E^ eading < 25 GeV the data lie above 

the predictions (see Fig. ||) and for events with very backward jets, — 1 < rf^ < 0, 
the measurement lies below the calculations (see Fig. ||). The Monte Carlo studies 
discussed in section |] show that fragmentation decreases the measured cross section for 
events with negative pseudorapidities. It is therefore to be expected that the NLO-QCD 
calculations, in which no parton-to-hadron fragmentation is included, predict a higher 
cross section than that observed in this region. 

The dijet cross section is also presented as a function of the pseudorapidity of one of the 
jets while keeping the other jet fixed in specific pseudorapidity ranges. Numerical values 
for the cross section and the uncertainties are given in Tab. [3] and || and are plotted in 
Fig. H The cross section peaks for events with rj^ et near 1 and falls rapidly for events 
with rf^ 1 < 0. The measurements are again compared to NLO-QCD calculations, but 
now using three different parameterisations for the parton densities in the photon. For 
the full x° bs range, at central and forward pseudorapidities of the jets, the data lie above 
all predictions. At backward pseudorapidities, as was the case for the cross section as a 
function of E^ eading , the data lie below the calculations. In the high x° bs region general 
agreement is seen between the data and the predictions. 

Fig. |5]d shows a comparison between the NLO-QCD results from four different groups 
for the range < r/j e< < 1. Each calculation uses the same parton density distributions 
for the proton, CTEQ4M p), and the photon, GRV-HO The calculations from 

Aurenche et al., Frixione et al., Harris et al. and Klasen et al. agree to within a few 
percent. 

In summary, it has been shown that NLO-QCD calculations generally describe the 
measured cross sections. However, for backward pseudorapidities the data are below 
the calculations, which is expected due to fragmentation, while for forward and central 
pseudorapidities the data are above the NLO predictions. In the latter kinematic region 
theoretical uncertainties are expected to be small. 



9.2 Cross sections for 212 < W 1P < 277 GeV 

The pseudorapidity dependence of the cross section has also been determined for events 
in a narrower region in y, which corresponds to a narrower range in W 7P , the photon- 
proton CM energy. In such a region the sensitivity to the photon structure is expected 
to be larger. This follows from the relation between y, x° bs and the pseudorapidities 
of the jets (see formula ^). Using a narrower range of y values implies that the cross 
section for specific pseudorapidities of the jets corresponds to a narrower range of x° bs 
values. It is natural to select a narrow region of high y values rather than a narrow 
region of low y values, since, in the latter case, events with low x° bs would fall out of 
the range of jet pseudorapidities, —1 < if et < 2. 

Using a range of 0.50 < y < 0.85, the cross section is presented as a function of the 
pseudorapidity of one of the jets while keeping the other jet fixed in a specific pseu- 
dorapidity bin. Values for the cross section and the uncertainties are given in Tab. || 
and |6] and are shown in Fig. |(| The cross section for this high-y region peaks at more 
backward pseudorapidities than the cross section for the full y range, as observed in 
a previous ZEUS study Q, and also the peak is more pronounced than for the full y 
range. This observation is consistent with the expected closer correlation between jf ets 
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and x° when the y range is restricted. The peak in the cross sections at backward 
pseudorapidities reflects the peak near x° bs ~ 1 in Fig. [l| and the tail towards positive 
pseudorapidities corresponds to low x° bs values. 

The measurements are again compared to NLO-QCD calculations using the GRV-HO, 
AFG-HO and GS96-HO parameterisations of the photon structure. The NLO predic- 
tions show an enhanced sensitivity to the choice of parameterisation for the photon 
structure. In particular in the region 1 < rj{ et < 2 there are clear differences in shape 
between the NLO predictions corresponding to different parton densities in the photon. 
In the most backward bins, where r^ 1 < —0.5 or where rf^ < 0, the data again lie below 
the calculations, but, as stated above, fragmentation effects are large in this region. At 
central and forward pseudorapidities, both for the full and for the high x^ s range, the 
data lie above the NLO calculations. 

In Fig. ^d a comparison is again made between the NLO-QCD results from different 
groups. The calculations agree to within a few percent. 

The fact that the cross sections, measured in the region where jets are produced at 
central and forward pseudorapidities and where theoretical uncertainties are expected 
to be small, lie above the NLO-QCD predictions, suggests that in this kinematic region 
the parton densities in the photon are too small in the available parameterisations. The 
disagreement between the data and the calculations is observed for the full x° bs range 
and to a lesser extent also for x° bs > 0.75. It is strongest at central pseudorapidities. 
This region corresponds to values of x 7 that lie roughly between 0.5 and 1. 



10 Summary and conclusions 

A measurement of dijet photoproduction, in the range 0.20 < y < 0.85, Q 2 < 1 GeV 2 , 
-1 < rf et < 2, E^ eadmg > 14 GeV and E^ second > 11 GeV, has been presented. Jets 
are defined in the hadronic final state by applying the /cy-clustering jet algorithm. The 
cross section has been compared to NLO-QCD predictions. 

For the full y region, 0.20 < y < 0.85, corresponding to 134 < W 1P < 277 GeV, the 
dijet cross section has been measured as a function of the transverse energy of the 
leading jet and as a function of the pseudorapidities of the jets. The dependence on 
the transverse energy of the leading jet is generally well described by the NLO-QCD 
calculations, although for events with two forward-going jets and Ej? leading < 25 GeV, 
the data lie above the NLO-QCD calculations. Also, the cross section as a function 
of the pseudorapidities of the jets lies above the NLO-QCD calculations at central and 
forward pseudorapidities. In the region x° bs > 0.75, the calculations agree with the 
measured cross section. 

In the high-y region, 0.50 < y < 0.85 (212 < W 1P < 277 GeV), where a stronger 
sensitivity to the photon structure is expected, the cross section at central and forward 
pseudorapidities lies further above the predictions than for the full y range. Also the 
cross section lies above the NLO-QCD calculations for x° bs > 0.75. 
Since theoretical uncertainties are expected to be small in most of the kinematic regime 
of the present analysis, as was discussed in sections ||] and ^, the discrepancies observed 
between the data and the NLO-QCD calculations suggest that, in the kinematic region 
of the present analysis, the available parameterisations of the parton densities in the 
photon are too small. 
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The results presented in this paper cover a kinematic region where both x° and £jf , 
which acts as the factorisation scale, are high. This region has not been studied in 
measurements. It remains to be established whether the parton density functions in the 
photon can be modified to describe the present data while remaining consistent with the 
existing i 7 ^ data from e + e _ experiments. It is hoped that phenomenologists carrying 
out comprehensive NLO-QCD fits will be able to include the data in this paper in their 
fits to determine the parton density functions in the photon and thereby clarify this 
issue. 
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Figure 1: The x° bs spectrum of the selected dijet sample, compared to the HERWIG 5.9 and the 
PYTHIA 5.7 Monte Carlo predictions, which have been weighted as described in section |^. The direct 
component from the HERWIG Monte Carlo is shown separately as the shaded histogram. Only statistical 
uncertainties are plotted 
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Figure 2: The transverse energy flow around the jet axis (integrated over |A0| < 1), in three ranges 
of the transverse energy of the jet and in four bins in x° . The data are compared to the HERWIG 5.9 
predictions. For the data only statistical uncertainties are shown. 
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Figure 3: Dijet cross section as a function of E^ eading for rj{ et between 1 and 2, in three regions of rff* . 
The results for — 1 < rj!^ 1 < and < rj^ 1 < 1 are scaled by the factors indicated in the figure. The filled 
circles correspond to the entire x° B range while the open circles correspond to events with x° bs > 0.75. 
The shaded band indicates the uncertainty related to the energy scale. The thick error bar indicates the 
statistical uncertainty and the thin error bar indicates the systematic and statistical uncertainties added 
in quadrature. The data are compared to NLO-QCD calculations, using the GRV-HO parameterisation 
for the photon structure. 
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Figure 4: Dijet cross section as a function of E^f eadi . For the two upper sets of data T] 3 ^ 1 lies 
between and 1 and for the lower set of data rrj et lies between -1 and 0. The r^ et regions are indicated 
the figure. The two lower sets of data are scaled by the factors indicated in the figure. The filled circles 
correspond to the entire x° bs range while the open circles correspond to events with x° bs > 0.75. The 
shaded band indicates the uncertainty related to the energy scale. The thick error bar indicates the 
statistical uncertainty and the thin error bar indicates the systematic and statistical uncertainties added 
in quadrature. The data are compared to NLO-QCD calculations, using the GRV-HO parameterisation 
for the photon structure. 
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NLO-QCD, AFG-HO, x;">0.75 
NLO-QCD, GS96-HO, xf>0.75 



Figure d) '. 

HARRIS et ol„ GRV-HO 

KLASEN et al., GRV-HO 

FRIXIONE et al., GRV-HO 

AURENCHE et al., GRV-HO 



Figure 5: Figures a), b) and c) show the dijet cross section as a function of r/^ et in bins of rj\ et . The filled 
circles correspond to the entire x° ba range while the open circles correspond to events with x° bB > 0.75. 
The shaded band indicates the uncertainty related to the energy scale. The thick error bar indicates 
the statistical uncertainty and the thin error bar indicates the systematic and statistical uncertainties 
added in quadrature. The full, dotted and dashed curves correspond to NLO-QCD calculations, using 
the GRV-HO, GS96-HO and the AFG-HO parameterisations for the photon structure, respectively. In 
d) the NLO-QCD results for the cross section when < < 1 and for a particular parameterisation 
of the photon structure are compared. 
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Figures a), b) and c); 

a ZEUS 1995 o ZEUS 1995, xf*>0. 75 

NLO-QCD, GRV-HO NLO-QCD, GRV-HO, xf>0.75 

NLO-QCD, AFG-HO NLO-QCD, AFG-HO, xf >0.75 

NLO-QCD, GS96-HO NLO-QCD, GS96-H0, xf* >0.75 



Figure d); 

HARRIS et al., GRV-HO 

KLASEN et al„ GRV-HO 

FRIXIONE et al., GRV-HO 



Figure 6: Figures a), b) and c) show the dijet cross section as a function of rf 2 et in bins of r^ et and for 
0.50 < y < 0.85. The filled circles correspond to the entire x° bs range while the open circles correspond 
to events with x° ba > 0.75. The shaded band indicates the uncertainty related to the energy scale. The 
thick error bar indicates the statistical uncertainty and the thin error bar indicates the systematic and 
statistical uncertainties added in quadrature. The full, dotted and dashed curves correspond to NLO- 
QCD calculations, using the GRV-HO, GS96-HO and the AFG-HO parameterisations for the photon 
structure, respectively. In d) the NLO-QCD results for the cross section when < rf^ 1 < 1 and for a 
particular parameterisation of the photon structure are compared. 
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da/dE J T et leadinq for: 0.20 < y < 0.85 and all if values 


E jet 

T leading 

GeV 


da / dE T leading 

pb/GeV 


^■stat 

pb/GeV 


&syst (+/-) 
pb/GeV 


A-E-scale {+ / ~) 

pb/GeV 






-1 < 7]{ et < & - 


1 < rf 2 et < 




14.0 .. 


17.0 


13.6 


1.2 


0.9/ -2.2 


3.2/ -1.2 


17.0 .. 


21.0 


1.81 


0.37 


1.16/ -0.54 


0.60/ -0.26 






u < Vi 


<1&-1< rf 2 et < 




14.0 .. 


17.0 


46.2 


1.7 


4.4/ -8.9 


4.1/ -2.1 


17.0 .. 


21.0 


18.3 


0.9 


0.7/ -3.1 


2.7/ -1.5 


21.0 .. 


25.0 


4.4 


0.4 


0.6/ -0.2 


0.8/ -0.5 


25.0 .. 


29.0 


0.97 


0.20 


0.29/ -0.29 


0.32/ -0.15 


29.0 .. 


35.0 


0.093 


0.046 


0.030/-0.063 


0.035/-0.023 






1 < rff 


< 2 & - 1< rf 2 et < 




14.0 .. 


17.0 


41.2 


1.6 


2.4/ -7.1 


5.8/ -1.0 


17.0 .. 


21.0 


16.9 


0.8 


0.5/ -1.4 


1.2/ -1.2 


21.0 .. 


25.0 


5.1 


0.5 


0.7/ -0.4 


0.7/ -0.6 


25.0 .. 


29.0 


1.56 


0.26 


0.17/ -0.23 


0.29/ -0.19 


29.0 .. 


35.0 


0.42 


0.11 


0.24/ -0.05 


0.10/ -0.08 


< r?f * < 1 k < rf> 2 et < 1 


14.0 .. 


17.0 


81.8 


3.0 


3.1/ -1.8 


11.9/ -1.6 


17.0 .. 


21.0 


42.5 


1.9 


0.5/ -4.8 


4.1/ -2.6 


21.0 .. 


25.0 


18.2 


1.2 


2.4/ -2.5 


2.6/ -1.1 


25.0 .. 


29.0 


7.5 


0.8 


0.3/ -0.7 


1.1/ -0.8 


29.0 .. 


35.0 


2.4 


0.4 


0.1/ -0.3 


0.3/ -0.3 


35.0 .. 


41.0 


0.49 


0.16 


0.14/ -0.04 


0.12/ -0.05 


1 < rfi * < 2 & < rf 2 et < 1 


14.0 .. 


17.0 


73.7 


2.0 


1.7/ -3.0 


9.0/ -0.5 


17.0 .. 


21.0 


40.4 


1.3 


0.7/ -2.6 


4.7/ -2.9 


21.0 .. 


25.0 


17.9 


0.9 


0.2/ -1.3 


1.9/ -1.6 


25.0 .. 


29.0 


8.2 


0.6 


0.1/ -1.0 


1.1/ -1.0 


29.0 .. 


35.0 


2.8 


0.3 


0.7/ -0.6 


0.4/ -0.4 


35.0 .. 


41.0 


1.18 


0.18 


0.25/ -0.41 


0.13/ -0.14 


41.0 .. 


48.0 


0.20 


0.07 


0.15/ -0.03 


0.04/ -0.02 


48.0 .. 


55.0 


0.28 


0.10 


0.04/ -0.19 


0.06/ -0.05 


1 < 7]{ et < 2 & 1< rf 2 1 < 2 


14.0 .. 


17.0 


49.6 


2.3 


4.0/ -2.0 


7.1/ -0.8 


17.0 .. 


21.0 


30.4 


1.6 


1.0/ -2.9 


3.4/ -2.7 


21.0 .. 


25.0 


15.0 


1.1 


0.5/ -1.9 


1.6/ -1.5 


25.0 .. 


29.0 


6.2 


0.7 


0.7/ -0.7 


0.8/ -0.6 


29.0 .. 


35.0 


2.8 


0.4 


0.3/ -0.3 


0.3/ -0.4 


35.0 .. 


41.0 


1.53 


0.29 


0.05/ -0.60 


0.25/ -0.12 


41.0 .. 


48.0 


0.39 


0.14 


0.05/ -0.06 


0.06/ -0.07 


48.0 .. 


55.0 


0.099 


0.070 


0.183/-0.009 


0.009/-0.007 



Table 1: The dijet cross section for the full x° ba range and 0.20 < y < 0.85, as a function of E^ eadin 
in bins of the jet pseudorapidities. 
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da/dE J T et leading for: 


0.20 < y < 0.85 and xf 


> U. 1 


E jet 

T leading 

GeV 


da / dE T leading 

pb/GeV 


^■stat 

pb/GeV 


&syst (+/-) 
pb/GeV 


^E-scale \+ / ~ ) 

pb/GeV 






-1 < 7]{ et < & - 


1 < rf 2 et < 




14.0 .. 


17.0 


12.4 


1.2 


1.0/ -2.0 


2.9/ -1.1 


17.0 .. 


21.0 


1.74 


0.36 


1.12/ -0.52 


0.58/ -0.25 






u < Vi 


<1&-1< rf 2 et < 




14.0 .. 


17.0 


37.1 


1.5 


3.5/ -7.7 


3.3/ -1.7 


17.0 .. 


21.0 


15.4 


0.8 


1.1/ -2.5 


2.3/ -1.3 


21.0 .. 


25.0 


4.2 


0.4 


0.5/ -0.2 


0.8/ -0.5 


25.0 .. 


29.0 


0.93 


0.19 


0.22/ -0.32 


0.30/ -0.14 


29.0 .. 


35.0 


0.093 


0.046 


0.030/-0.063 


0.035/-0.023 






1 < rff 


< 2 & - 1< rf 2 et < 




14.0 .. 


17.0 


26.2 


1.2 


3.2/ -5.9 


3.7/ -0.7 


17.0 .. 


21.0 


11.8 


0.7 


0.7/ -1.1 


0.9/ -0.9 


21.0 .. 


25.0 


4.1 


0.4 


0.6/ -0.5 


0.6/ -0.5 


25.0 .. 


29.0 


1.48 


0.25 


0.13/ -0.22 


0.27/ -0.18 


29.0 .. 


35.0 


0.39 


0.11 


0.20/ -0.02 


0.10/ -0.07 


< r?f * < 1 k < rf> 2 et < 1 


14.0 .. 


17.0 


48.6 


2.3 


4.4/ -2.1 


7.1/ -0.9 


17.0 .. 


21.0 


27.8 


1.5 


1.6/ -3.5 


2.7/ -1.7 


21.0 .. 


25.0 


13.2 


1.0 


2.1/ -2.0 


1.9/ -0.8 


25.0 .. 


29.0 


6.1 


0.7 


0.1/ -1.0 


0.9/ -0.7 


29.0 .. 


35.0 


1.9 


0.3 


0.2/ -0.3 


0.2/ -0.3 


35.0 .. 


41.0 


0.49 


0.16 


0.11/ -0.10 


0.12/ -0.05 


1 < rfi * < 2 & < rf 2 et < 1 


14.0 .. 


17.0 


29.0 


1.3 


1.1/ -3.4 


3.6/ -0.2 


17.0 .. 


21.0 


18.4 


0.9 


1.6/ -1.9 


2.1/ -1.3 


21.0 .. 


25.0 


8.7 


0.6 


0.6/ -0.7 


0.9/ -0.8 


25.0 .. 


29.0 


5.1 


0.4 


0.1/ -0.8 


0.7/ -0.6 


29.0 .. 


35.0 


1.85 


0.23 


0.39/ -0.35 


0.30/ -0.26 


35.0 .. 


41.0 


0.83 


0.15 


0.30/ -0.32 


0.09/ -0.10 


41.0 .. 


48.0 


0.125 


0.056 


0.092/-0.012 


0.023/-0.012 


48.0 .. 


55.0 


0.21 


0.09 


0.06/ -0.16 


0.05/ -0.04 


1 < 7]{ et < 2 & 1< rf 2 1 < 2 


14.0 .. 


17.0 


1.28 


0.37 


1.32/ -0.73 


0.18/ -0.02 


17.0 .. 


21.0 


4.2 


0.6 


0.6/ -1.2 


0.5/ -0.4 


21.0 .. 


25.0 


4.1 


0.6 


0.3/ -0.5 


0.4/ -0.4 


25.0 .. 


29.0 


2.8 


0.5 


0.3/ -0.5 


0.4/ -0.3 


29.0 .. 


35.0 


1.29 


0.28 


0.38/ -0.21 


0.16/ -0.16 


35.0 .. 


41.0 


0.91 


0.23 


0.05/ -0.32 


0.15/ -0.07 


41.0 .. 


48.0 


0.24 


0.11 


0.07/ -0.06 


0.04/ -0.04 


48.0 .. 


55.0 


0.099 


0.070 


0.096/-0.009 


0.009/-0.007 



Table 2: The dijet cross section for x° bs > 0.75 and 0.20 < y < 0.85, as a function of E^\ eaiing in bins 
of the jet pseudorapidities. 
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da/d^ et for: 0.20 < y < 0.85 and all if values 


■let 

m 


da/dif 2 et 




&syst (+/") 


A-E-scale (+/ — ) 






pb 


pb 


pb 


pb 








-1 < 


rff < 




-0.5 .. 


0.0 


88 


5 


6/ -13 


20/ -9 


0.0 .. 


0.5 


209 


9 


10/ -39 


26/ -13 


0.5 .. 


1.0 


258 


9 


4/ -34 


28/ -17 


1.0 .. 


1.5 


240 


9 


9/ -39 


28/ -12 


1.5 .. 


2.0 


201 


8 


7/ -17 


25/ -10 


< r?f * < 1 


-1.0 .. 


-0.5 


115 


7 


4/ -30 


21/ -12 


-0.5 .. 


0.0 


353 


11 


11/ -44 


35/ -19 


0.0 .. 


0.5 


513 


13 


16/ -29 


66/ -25 


0.5 .. 


1.0 


558 


14 


19/ -43 


71/ -23 


1.0 .. 


1.5 


541 


13 


1/ -15 


61/ -32 


1.5 .. 


2.0 


486 


13 


6/ -23 


63/ -25 


1 < 7]{ et < 2 


-1.0 .. 


-0.5 


113 


6 


10/ -15 


18/ -8 


-0.5 .. 


0.0 


328 


11 


7/ -42 


36/ -15 


0.0 .. 


0.5 


479 


12 


4/ -19 


60/ -25 


0.5 .. 


1.0 


549 


14 


1/ -18 


63/ -33 


1.0 .. 


1.5 


416 


12 


9/ -14 


54/ -28 


1.5 .. 


2.0 


358 


11 


12/ -15 


43/ -23 



Table 3: The dijet cross section, for all x° ba values and 0.20 < y < 0.85, as a function of rf 2 et , for rf x 
fixed. 
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da/drf et for: 0.20 < y 


< 0.85 and x° bs > 0.75 


jet 

m 




da/dif 2 et 




&syst (+/") 


A-E-scale (+/ — ) 






pb 


pb 


pb 


pb 








-1 < 


rff < 




-0.5 .. 


0.0 


80 


5 


8/ -12 


18/ -8 


0.0 .. 


0.5 


185 


8 


13/ -35 


23/ -11 


0.5 .. 


1.0 


204 


8 


9/ -30 


23/ -13 


1.0 .. 


1.5 


173 


8 


15/ -32 


20/ -9 


1.5 .. 


2.0 


129 


7 


12/ -19 


16/ -7 


< r?f * < 1 


-1.0 .. 


-0.5 


109 


7 


3/ -29 


20/ -11 


-0.5 .. 


0.0 


283 


10 


14/ -37 


28/ -15 


0.0 .. 


0.5 


359 


11 


19/ -25 


46/ -18 


0.5 .. 


1.0 


339 


11 


26/ -35 


43/ -14 


1.0 .. 


1.5 


273 


9 


15/ -27 


31/ -16 


1.5 .. 


2.0 


195 


8 


14/ -20 


25/ -10 


1 < 7]{ et < 2 


-1.0 .. 


-0.5 


94 


6 


11/ -13 


15/ -7 


-0.5 .. 


0.0 


210 


8 


14/ -35 


23/ -9 


0.0 .. 


0.5 


241 


9 


20/ -27 


30/ -12 


0.5 .. 


1.0 


227 


9 


7/ -18 


26/ -14 


1.0 .. 


1.5 


95 


6 


6/ -8 


12/ -6 


1.5 .. 


2.0 


30 


3 


2/ -2 


4/ -2 



Table 4: The dijet cross section for x°f s > 0.75 and 0.20 < y < 0.85, as a function of rf 2 et , for rft 
fixed. 
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da/d^ et for: 0.50 < y < 0.85 and all if values 


■let 

m 


da/drf 2 et 




&syst (+/") 


A-E-scale (+/ — ) 






pb 


pb 


pb 


pb 








-1 < 


rff < 




-0.5 .. 


0.0 


88 


5 


6/ -13 


20/ -9 


0.0 .. 


0.5 


208 


9 


12/ -41 


25/ -11 


0.5 .. 


1.0 


232 


9 


2/ -35 


25/ -14 


1.0 .. 


1.5 


185 


8 


11/ -26 


20/ -8 


1.5 .. 


2.0 


152 


7 


2/ -19 


14/ -8 


< r?f * < 1 


-1.0 .. 


-0.5 


115 


7 


4/ -30 


21/ -12 


-0.5 .. 


0.0 


326 


11 


11/ -47 


29/ -15 


0.0 .. 


0.5 


284 


10 


6/ -13 


48/ -14 


0.5 .. 


1.0 


218 


9 


3/ -13 


32/ -13 


1.0 .. 


1.5 


162 


7 


12/ -2 


26/ -10 


1.5 .. 


2.0 


153 


7 


3/ -3 


21/ -8 


1 < 7]{ et < 2 


-1.0 .. 


-0.5 


110 


6 


10/ -15 


16/ -8 


-0.5 .. 


0.0 


227 


9 


8/ -32 


22/ -14 


0.0 .. 


0.5 


186 


8 


4/ -4 


29/ -9 


0.5 .. 


1.0 


128 


6 


11/ -3 


18/ -8 


1.0 .. 


1.5 


122 


6 


4/ -7 


17/ -8 


1.5 .. 


2.0 


117 


6 


11/ -5 


17/ -5 



Table 5: The dijet cross section, for all x° ba values and 0.50 < y < 0.85, as a function of rf 2 et , for rf x 
fixed. 
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dcr/dri iet for: 0.50 < y 


< 0.85 and x° bs > 0.75 


let 

m 




da/drf 2 et 


^■stat 


A-syst (+/") 


&-E-scale (+/ — ) 






pb 


pb 


pb 


pb 


-1 < 7]{ £t < 


-0.5 .. 


0.0 


80 


5 


8/ -12 


18/ -8 


0.0 .. 


0.5 


183 


8 


15/ -37 


22/ -10 


0.5 .. 


1.0 


178 


8 


4/ -29 


19/ -11 


1.0 .. 


1.5 


122 


6 


13/ -18 


13/ -5 


1.5 .. 


2.0 


84 


5 


3/ -15 


8/ -4 


< 7]{ et < 1 


-1.0 .. 


-0.5 


108 


7 


4/ -29 


20/ -11 


-0.5 .. 


0.0 


255 


9 


9/ -38 


23/ -12 


0.0 .. 


0.5 


152 


7 


5/ -9 


26/ -8 


0.5 .. 


1.0 


66 


5 


0/ -7 


10/ -4 


1.0 .. 


1.5 


27 


3 


2/-1 


4/ -2 


1.5 .. 


2.0 


15 


2 


2/ -1 


2/-1 


1 < r?f < 2 


-1.0 .. 


-0.5 


91 


6 


10/ -13 


13/ -6 


-0.5 .. 


0.0 


118 


6 


5/ -18 


12/ -7 


0.0 .. 


0.5 


37 


3 


0/ -4 


6/ -2 


0.5 .. 


1.0 


5.0 


1.2 


2.5/ -0.5 


0.7/ -0.3 



Table 6: The dijet cross section for x° b3 > 0.75 and 0.50 < y < 0.85, as a function of rf 2 et , for r){ 
fixed. 
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